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Abstract 

Self-avoiding walks self-interacting via nearest neighbours (ISAW) and self-avoiding trails in- 
teracting via multiply-visited sites (ISAT) are two models of the polymer collapse transition of 
a polymer in dilute solution. On the square lattice it has been established numerically that the 
collapse transition of each model lies in a different universality class. 

It has been shown that by adding stiffness to the ISAW model a second low temperature phase 
eventuates and a more complicated phase diagram ensues with three types of transition that meet 
at a multi-critical point. For large enough stiffness the collapse transition becomes first-order. 
Interestingly, a phase diagram of a similar structure has been seen to occur in an extended ISAT 
model on the triangular lattice without stiffness. It is therefore of interest to see the effect of adding 
stiffness to the ISAT model. 

We have studied by computer simulation a generalised model of self-interacting self-avoiding 
trails on the square lattice with a stiffness parameter added. Intriguingly, we find that stiffness 
does not change the order of the collapse transition for ISAT on the square lattice for a very wide 
range of stiffness weights. While at the lengths considered there are clear bimodal distributions 
for very large stiffness, our numerical evidence strongly suggests that these are simply finite-size 
effects associated with a crossover to a first-order phase transition at infinite stiffness. 



1 Introduction 



The collapse transition of a polymer in a dilute solution has been a continuing focus of study in lattice 
statistical mechanics for decades pi, f2]. This transition describes the change in the scaling of the 
polymer with length that occurs as the temperature is lowered. At high temperatures the radius of 
gyration of polymer scales in a way swollen relative to a random walk: this is known as the excluded 
volume effect. At low temperatures a polymer condenses into dense, usually disordered, globule, with 
a much smaller radius of gyration. The interest in this phase transition has occurred both because of 
the motivation of physical systems but also because of the study of integrable cases [3l H] of lattice 
models, that have proved especially fruitful in two dimensions. While the canonical lattice model of 
the configurations of a polymer in solution has been the self-avoiding walk (SAW), where a random 
walk on a lattice is not allowed to visit a lattice site more than once, an alternative has been to use 
bond-avoiding walks, or a self-avoiding trail. A self-avoiding trail (SAT) is a lattice walk configuration 
where the excluded volume is obtained by preventing the walk from visiting the same bond, rather 
than the same site, more than once. These were used initially to model polymers with loops [5j but 
have subsequently occurred in integrable loop models in two dimensions [Ij. A model of collapsing 
polymers can be constructed starting from self-avoiding trails, known as interacting self-avoiding trails 
(IS AT). Here energies are associated with multiply- visited sites and by favouring configurations with 
many such sites a collapse transition can be initiated. 

Owczarek and Prellberg studied numerically the ISAT collapse on the square lattice by two different 
approaches [6l [7] and in either case found a strong continuous transition with specific heat exponent 
a = 0.81(3). Recently, on the triangular lattice Doukas et al. [8] found that by changing the weighting 
of doubly and triply visited sites a first-order transition can ensure or alternatively, depending on 
the ratio of these weightings, a weaker second-order transition that mimics the collapse found in the 
canonical interacting self-avoiding walk (ISAW) model (also know as the ^-point). They also found 
that the low temperature phase could become fully dense rather than globular. 

Coorespondingly, there is also a modification of the ISAW model that displays two phase transi- 
tions for a range of parameters, namely the semi-flexible ISAW model [9t 110^ [TT]. Here two energies 
are included: the nearest-neighbour site interaction of the ISAW model and also a stiffness energy 
associated with consecutive parallel bonds of the walk (equivalently, a bending energy for bends in 
the walk). This has been studied on the cubic lattice by Bastolla and Grassberger |9j. They showed 
that when there is a strong energetic preference for straight segments, this model undergoes a single 
first-order transition from the excluded-volume high-temperature state to a fully dense state. On the 
other hand, if there is only a weak preference for straight segments, the polymer undergoes two phase 
transitions. On lowering the temperature the polymer undergoes a 0-point transition to the liquid 
globule followed by a first-order transition to the fully dense phase at a lower temperature. Recent 
work by Krawczyk et al. [12] concerning the ISAW model on the square lattice in presence of a stiffness 
parameter showed that the introduction of stiffness can change the universality class of the collapse 
transition in two dimensions. For large stiffness the transition becomes first order and the collapsed 
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phase moves from being globular to fully dense. 

Recently, Foster [13] introduced and studied a generalised ISAT model on the square lattice which 
incorporates stiffness. Using transfer matrices and the phenomenological renormalisation group, that 
study predicted that the ISAT universality class is unaffected by a range of values of stiffness. However, 
the results suggested the appearance of a first-order transition for sufficiently large stiffness. 

In this work we use Monte Carlo simulation to explore ISAT in presence of stiffness and the 
predictions of Foster [13] . We also explore the low temperature phase of the model and find that there 
is only one low temperature phase and that it is fully dense for the range of stiffness studied. 



2 ISAT 

The model of interacting trails on the square lattice is defined as follows. Consider the ensemble Tn 
of self-avoiding trails (SAT) of length n, that is, of all lattice paths of n steps that can be formed on 
the square lattice such that they never visit the same bond more than once. Given a SAT ipn ^ Tn, 
we associate an energy —et with each doubly visited site, and denote their number by m{tpn)- The 
probability of il^n is then given by 

p/3etm(i/)„) 

(2.1) 



zlsAT^ry 

where we define the Boltzmann weight uot = exp(/3et) and /3 is the inverse temperature l/ksT. The 
partition function of the ISAT model is given by 



Zf^^(r)= ^ uj'^^^"\ (2.2) 



The finite-length reduced free energy is 

K„(r) = -log Z„(T) (2.3) 
n 

and the thermodynamic limit is obtained by taking the limit of large n, i.e., 

k{T) = lim Kn{T). (2.4) 

n— ^-oo 

It is expected that there is a collapse phase transition at a temperature Tc characterized by a non- 
analyticity in k{T). 

The collapse transition can be characterized via a change in the scaling of the size of the polymer 
with temperature. It is expected that some measure of the size, such as the radius of gyration or the 
mean squared distance of a monomer from the end points, i?^(T), scales at fixed temperature as 

Rl{T) ~ An^" (2.5) 

with some exponent ly. At high temperatures the polymer is swollen and in two dimensions it is 
accepted that = 3/4 [[3J. At low temperatures the polymer becomes dense in space, though not 
necessarily space filling, and the exponent is = 1/2. However, for the ISAT model the collapsed 
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phase has been seen to be space fiUing [14J. If the cohapse transition is second-order, the scahng at 
Tc of the size is intermediate between the high and low temperature forms. In the thermodynamic 
hmit the expected singularity in the free energy can be seen in its second derivative (the specific heat). 
Denoting the (intensive) finite length specific heat per monomer by Cn{T), the thermodynamic limit 
is given by the long length limit as 

C{T) = lim CniT) . (2.6) 

n— >oo 

One expects that the singular part of the specific heat behaves as 

C(T) ~ BlTc -T|-" , (2.7) 

where a < 1 for a second-order phase transition. The singular part of the thermodynamic limit 
internal energy behaves as 

[/(T) ~B|rc-r|^-" , (2.8) 

if the transition is second-order, and there is a jump in the internal energy if the transition is first-order 
(an effective value of a = 1). 

Moreover, one expects crossover scaling forms [T5] to apply around this temperature, so that 

c„(r)~n°'^c((r-r,)n*) (2.9) 

with < (/) < 1 if the transition is second-order, and 

Cn{T)^nC{{T-T^)n) (2.10) 

if the transition is first-order. From |15| we point out that the exponents a and (f) are related via 

1 

2-a = -. (2.11) 



Important for numerical estimation is the use of equation (2.9) at the peak value of the specific heat 
given by yP"'''' = (T - Tc)n'^ so that 

^^eakf^j.^ ^ (.peak ^ac^ (2.12) 

where C^^"*^ = C{y'P^°'^) is a constant. 

A previous study [7J of ISAT model on the square lattice has shown that there is a collapse 
transition with a strongly divergent specific heat, with 

Q(/) = 0.68(5) (2.13) 

and so the individual exponents have been estimated as 

(j) = 0.84(3) and a = 0.81(3) . (2.14) 

At T = Tc, given bv [TBI [TT] 

ut = 3, (2.15) 

it was predicted [6] that 

Rl{T)^An{\ognf . (2.16) 
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3 Semi-flexible ISAT 
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Figure 1: An example of semi- flexible ISAT configuration with three (m = 3) multiple-visit interactions 
associated with the Boltzmann weight r and seven (s = 7) straight segments each associated with the 
Boltzmann weight p. Note that when the trail crosses itself, the crossing site is associated with a total 
weight T]? as a crossing necessarily requires two straight segments. 



The semi-flexible ISAT (SFISAT) model can be defined as follows. Consider the set of bond- 
avoiding paths Tn as defined in the previous section. Given a SAT i^n ^ Tn, '^^ associate an energy 
—£t every time the path visits the same site more than once, as in ISAT. Additionally, we define a 
straight segment of the trail by two consecutive parallel edges, and we associate an energy —Eg to each 
straight segment of trail. 

For each configuration ipn ^ %i count the number m{ipn) of doubly- visited sites and s{ipn) 
of straight segments: see Figure [TJ Hence we associate with each configuration a Boltzmann weight 
^m{tpn) psi'tpn) ^i^Qj-Q J- — exp(/3et), p = exp(/3es), and /? is the inverse temperature l/ksT. The partition 
function of the SFISAT model is given by 



The probability of a configuration ipn is then 



p(V'n;T,p) 



(3.1) 



(3.2) 



Zn{T,p) 

When we set p = 1 the trail is fully flexible and the model reduces to the ISAT model. On the other 
end, if we set p = straight segments are excluded and our model requires the path to turn at every 
site: this is known as the "L-lattice". Trails on the L-lattice may be mapped [18] into the Interacting 
Self- Avoiding Walk model on the Manhattan lattice ^19j. As a consequence of the mapping, the 
transition of ISAT on the L-lattice has been shown to be 6-like with a convergent specific heat, unlike 
ISAT on the square lattice. 

The average of any quantity Q over the ensemble set of path Tn is given generically by 



{Q){n;T,p) 



Qilpn)pitpn;r,p). 



(3.3) 



In particular, we can define the average number of doubly-visited sites per site and their respective 
fluctuations as 

(m) (m^) — (m)^ 



u 



n 
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n 



(3.4) 



One can also consider the average number of straight sections of the trail and their fluctuations 

„2\ /„\2 



nW = M and (3.5) 
n n 

An important quantity for what follows is the proportion of the sites on the trail that are at lattice 
sites which are not doubly occupied: 

.. = 1-"^. (3.6) 

n 

Foster [TJ] predicted that the universality class of fully flexible ISAT at p = 1 extends to other values 
of p and also that for large p there may be a change to a first-order transition. 
Defining the collapse as occurring at Tc{p) for constant p, it follows that 

Tail) = 3 and rc(0) = 2 (3.7) 

from the ISAT and L-lattice [IHl US HZ] results. 

4 Results 

We began by simulating the full two parameter space by using the flatPERM algorithm [20J. Flat- 
PERM outputs an estimate Wn^k of the total weight of the walks of length n at fixed values of some 
vector of quantities k = (ki, k2, ■ ■ ■ ,k£). From the total weight one can access physical quantities over 
a broad range of temperatures through a simple weighted average, e.g. 



{0)n{r) = A4^^ • (4.1) 



The quantities kj may be any subset of the physical parameters of the model. In our case we begin 
by using ki = m and k2 = s. 

We have simulated SFISAT using the full two-parameter flatPERM algorithm up to length n = 512, 
with 10^ iterations, collecting 7.8 • 10^ samples at the maximum length. 

To obtain a landscape of possible phase transitions we plot the largest eigenvalue of the matrix of 
second derivatives of the free energy with respect r and p (measuring the fluctuations and covariance 
in m and s) at length n = 512 in Figure [2] 

We notice that the strong peak corresponding to the ISAT transition at {t,p) = (3,1) extends 
upward to larger values of p, becoming stronger as p increases, as predicted in |13j . At smaller values 
oi p the peak of the specific heat seems to get weaker as it reaches the point (r, ) = (2,0). At this 
point, corresponding to a 9-like transition, the specific heat is known to not diverge. 

To investigate the nature of the transition at larger values of the stiffness parameter p, we have 
run extensive simulations of the model at few fixed values of p all up to length n = 1024. We have 
simulated different values of p between p = 0.1 and p = 10, running between S ~ 2.3-10^ and S ~ 4-10^ 
iterations, and collecting between 6.0 • 10^ and 1.7 • 10^*^ samples at the maximum length. Following 
|20] . we also measured the number of samples adjusted by the number of their independent growth 
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Figure 2: Density plot of the logarithm of the largest eigenvalue Xmax of the matrix of second deriva- 
tives of the free energy with respect to r and p at length 512. Darker regions denote larger values. 

steps between S"^-^-^ ~ 4.6 • lO''' and 4.4 • 10^ "effective samples" at the maximum length. Additionally, 
for the sake of comparison, we simulated the model at p = by putting the trails on the L-lattice, with 
S = 10^ iterations we collected 2.2 • 10^ samples at the maximum length, corresponding to 7.4 • 10'' 
effective samples. 



4.1 Specific heat 

We have begun by analysing the scaling of the specific heat by calculating the location of its peak 
T* = argmax^ c„(r) and thereby evaluating c* = c„(r*). In Figure [sj we plot the peak values of the 
specific heat for some of the models we have simulated. The exponent associated with the peak of 



the specific heat, see equation (2.12), is a0 if the transition is second order. For p = 1 we estimate 
a(p ~ 0.64 from trail lengths up to 1024, which is a little less that our previous estimate of 0.68(5) 
based upon much longer length trails |7|. For p > 1 we estimate a(p ~ 0.74 at p = 5 and a(j) ~ 0.67 at 
p = 10, both of which are compatible with that previous estimate for the ISAT transition of 0.68(5). 
For p = we expect that the specific heat does not diverge, but at the lengths we consider it is not 
surprising to measure a weakly increasing specific heat peak with significant curvature in a log-log 
plot. We can then compare this to the situation when p = 0.5. Although our estimate of a(j) ~ 0.52 is 
somewhat smaller at p = 0.5 than at p = 1, it is still significantly larger than at p = and unless one 
is willing to postulate continuously changing exponents this suggests that the entire line up to (but 
excluding) p = is in the ISAT universality class. There could, of course, be curvature in the plot 
that we cannot observe at this length scale. However, we reinforce our conclusion by examining the 
low temperature phase below. 
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Figure 3: Plot of the logarithm of the peak value of the specific heat versus the logarithm of the length 
of the trails for p = (left) and p = 0.5,5,10 (right). On the right, straight lines fit the data well, 
whereas there is clear curvature in the data on the left. 

Before we consider the low temperature behaviour of our model, we also observe that the location 
of the transition in the thermodynamic limit, lim„_i.oo T*{p) = Tc{p), seems to obey 

1 < rdp) < 3 for ah p, (4.2) 

attaining the value of Tc{p) = 3 when p = 1. As p increases from zero, Tc{p) increases until p reaches 
some value near one, and then decreases again. 

As indicated in Figure [4j for fixed large values of p, the location of the transition Tc{p) is close to, 
but larger than one. We therefore conjecture that Tc{p) — )• 1 as p — )• oo. The limit can be understood 
by the following ground state- type argument. For very large p the ground state depends on whether 
T is smaller or larger than one. For r < 1 the ground state is a single straight rod, while for r > 1 
it is a configuration that fills every edge of the lattice in order to maximise doubly visited sites (that 
is, a set of crossing long rods that join up on the surface of the polymer). We point out that such an 
argument also naturally leads to a first-order transition at infinite p. 



4.2 Energy distribution 

We then investigated the claim |13j that for large finite values of p the collapse transition becomes 
first-order. We did so by looking for the presence of bimodality in the energy distribution. In Figure [5] 
we plot the energy distribution in a neighbourhood of the critical point for p = 5 and p = 10. We 
don't find any evidence of a bimodality in either case. 

On the other hand, we do find evidence of a double peak forming at very large values of p. In 
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^0%00 0.001 0.002 0.003 0.004 

Figure 4: Plot of the location of the peak of the specific heat for p = 100 versus inverse length 1/n. 
A simple linear extrapolation gives us the estimate Tc( 100) ~ 1.1. 

Figure [6] we plot the energy distribution near the transition for p = 100, which has a clear double 
peak. Given the lengths of trails considered it is not clear immediately whether this is a manifestation 
of a crossover to an infinite p behaviour, which we argued above may well be first-order, or a real 
change in the order of the collapse transition at finite p. 

To investigate further we considered the smallest value of p at which there is a non-zero latent heat 
at the transition point T*{p) as a function of n: we denote this as pi{n). In Figure [7] we plot pi(n). 

We immediately see that pi{n) increases in n. Hence, if the transition becomes first-order at 
some finite value of p for infinite length trails, this would only occur for very large values of p. In 
fact, one can see that the increase in pi{n) is compatible with a power law such as n^/^, which one 
would expect if the finite-size corrections were related to surface contributions. Additionally, we have 
considered the finite-size latent heat at fixed large p. An extrapolation against n is compatible with 
a thermodynamic value of zero. This scenario would imply that we are simply seeing the crossover 
to the infinite-p behaviour in the simulations, and that the transition in the thermodynamic limit of 
infinite length stays second-order for all finite values of p. 

4.3 Low temperature region 

There is strong evidence that the low-temperature phase of the ISAW model is a globular phase that 
is not fully dense, while for interacting trails the low-temperature phase is maximally dense [U], i.e., 
the trail fills the lattice asymptotically. Therefore, for trails the portion of steps not involved with 
doubly-visited sites should tend to zero as n — )• oo in this phase. Following the analysis in [51 [T3] 
we measured the proportion Vn of steps not visiting the same site twice at r = 5, corresponding 
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Figure 5: Energy distribution across the critical region at length n = 1024 for p = 5 (a) and p = 10 
(b). For both figures the chosen temperatures are rieft = 1.85, Tmiddie = 1-93, and Tright = 2.02. The 
weight Tjniddie Corresponds to the location of the specific heat peak at this length (r* in the text), 
while Tieft and Tright are given by Tmiddie i ^At^ where At„ is the width of the specific heat peak at 
that length as measured at its half height. The distributions are normalised to have equal area. 

to a sufficiently low temperature to be in the collapsed phase. Figure [8] shows a plot of Vn against 
n~^/^, clearly indicating that f „ — )• within error bars. This reinforces our earlier conclusion that 
for < p < 1 the nature of collapse is that of the standard fully-flexible ISAT model, rather than 
changing to the ISAW-like 0-transition that occurs when p = 0. 

5 Phase diagram and Conclusions 

We have studied a generalised model of semi-flexible interacting trails (SFISAT) by including a stiffness 
parameter. From our analysis, we conjecture that the ISAT universality class is unaffected by the 
presence of stiffness. The universality class only changes in the singular limits of p — )• and p ^ oo. 
In the former limit the transition is 0-like, whereas the transition turns flrst-order in the latter. 

While at the lengths considered there are clear bimodal distributions for large values of p, our 
numerical evidence strongly suggests that these are likely to be flnite-size effects associated with a 
crossover to a flrst-order phase transition at inflnite stiffness. 

Our results indicate that for all finite values of p the low-temperature phase is maximally dense as 
in the fully-flexible ISAT model. Hence, unlike the ISAW model, there continues to be only a single 
low temperature phase when stiffness is added to ISAT. Putting the conclusions together suggests the 
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0.0 0.1 0.2 0.3 0.4 0.5 



m/n 

Figure 6: Energy distribution at length n = 512 for p = 100 at Boltzmann weight r = 2.01, chosen 
so that the heights of the two peaks are roughly equal, very close to the specific heat peak located at 
T* = 1.98. 




Figure 7: Plot of smallest p value pi (n) for which an observable latent heat is encountered for a finite 
length trail, versus ^/n. This scale is chosen to reflect the possibility of finite-size corrections due to 
the presence of a surface in the bulk. 

phase diagram shown in Figure [9j 

It would be of some interest to examine the effect of stiffness on the canonical ISAT model on the 
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Figure 8: Plots of Vn, the proportion of steps visiting the same site once at r = 5, which is in the low 
temperature phase, against n~^/^. The scale chosen is the natural low temperature scale. 




Figure 9: Schematic diagram of the SFISAT parameter space. The filled circle (blue online) on the 
horizontal axes depicts the location of the 0-like transition at (r,p) = (2,0), and the filled square 
(red online) depicts the location of the fully-flexible ISAT transition at (t, p) = (3, 1). The solid line 
corresponds to an ISAT-like transition. The whole region to the right of the curve, excluding the p = 
axis, is maximally-dense. 



triangular lattice where the low temperature phase is globular rather than fully dense as it is on the 
square lattice considered in this paper. 
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